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Abstract 

At a future linear collider very precise measurements, typically with errors 
of < 1%, are expected to be achievable. Such an accuracy gives sensitivity 
to the quantum corrections, which therefore must be incorporated in theo- 
retical calculations in order to determine the underlying new physics param- 
eters from prospective linear collider measurements. In the context of the 
chargino-neutralino sector of the minimal supersymmetric standard model, 
this involves fitting one-loop predictions to prospective measurements of the 
cross sections, forward-backward asymmetries and of the accessible chargino 
and neutralino masses. Taking recent results from LHC SUSY and Higgs 
searches into account we consider three benchmark scenarios, each with char- 
acteristic features. Our analysis shows how an accurate determination of the 
desired parameters is possible, providing in addition access to the stop masses 
and mixing angle. 
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1 Introduction 



A linear collider (LC) [T}|3] will be an ideal environment for high preci- 
sion studies of physics beyond the standard model (BSM). A particularly 
well-motivated BSM theory is the Minimal Supersymmetric Standard Model 
(MSSM). This provides the lightest neutralino as a candidate to explain the 
evidence for WIMP-like dark matter in the universe |4j[5]. Further, natu- 
ralness arguments (see e.g. ref. [6]) support light higgsino-like charginos and 
neutralinos, as also predicted by GUT motivated SUSY models [7]. The 
charginos and neutralinos could therefore be within reach of a first stage 
linear collider. 

In order to precisely determine the nature of the underlying SUSY model, 
the fundamental parameters have to be determined in the most model- 
independent way possible. The determination of the U(l) parameter Mi, 
the SU(2) parameter M 2 , the higgsino parameter fi and tan/3, the ratio of 
the vacuum expectation values of the two neutral Higgs doublet fields, at the 
percent level via chargino and neutralino pair-production has been shown 
to be possible at LO (see ref. {§] and references therein). Due to the ex- 
pected high precision of mass and coupling measurements at the LC, as well 
as the fact that one-loop effects in the MSSM are often large, higher order 
effects have to be considered. Taking these corrections into account addi- 
tional MSSM parameters become relevant, such as the masses of the stops 
and sleptons, which are also so far weakly constrained by the LHC. 

In this paper we show how it would be possible to determine the fun- 
damental parameters of the chargino and neutralino sector at the LC, in- 
cluding the uncertainty arising due to higher order effects. Specifically, we 
calculate the next-to-leading order (NLO) corrections to the cross-section (er) 
and forward-backward asymmetry (Afb) for chargino production, and also 
to the chargino and neutralino masses. A number of next-to-leading order 
(NLO) calculations of chargino and neutralino pair production at the LC 
can be found in the literature [9-12 . We perform our calculations in the 



on-shell (OS) scheme such that, as far as possible, the physical interpreta- 
tion of the masses is clear. Recent work on the OS renormalization of the 



chargino- neutralino sector can be found in refs. 13 16 



By fitting experimental results to these loop corrected predictions we show 
that it is possible to extract the fundamental parameters of the MSSM La- 
gangian. However due to the greater number of parameters, performing the 
fit is more involved as for the LO analysis. Choosing three potential MSSM 
scenarios, we assess the impact of the loop corrections and the feasibility of 
such an extraction in each. We further investigate the impact of obtaining 
masses of the charginos and neutralinos from threshold scans rather than the 



1 



continuum (see ref. [i] ) on the resulting accuracy of the parameters obtained 
from the fit. 

The paper is organised as follows. In sec. [2] we introduce the process 
studied and define necessary notation. We then provide details of the calcu- 
lation of the loop corrections in sec. |3j including details of the renormalization 
scheme used. In sec. [4] we further discuss the method employed in order to 
fit to the MSSM parameters, define the scenarios considered, and present our 
results. Finally in sec. [5] we discuss the implications of the results of the fits. 



2 Process studied and tree-level relations 



In this paper we study the determination of the fundamental parameters 
in the chargino-neutralino sector of the MSSM, via chargino production at 
a LC. The charginos, x^, and neutralinos, x°, are the mass eigenstates of 
the gauginos and higgsinos, as seen from the relevant part of the MSSM 



Lagrangian 17 



P L (U*XV% - PniVX^U 1 



+ 5« 



P L (N*YN% - P R (NY^N T ) l3 )xl 



where Pl/r = 1/2(1 =F 7s)- The mass matrix for the charginos is given by 



X 



M 2 V2M wS/5 

y/2M W C/3 /! 



(2) 



where sp/cp = sin/3/ cos/3, and M w is the mass of the W boson. This matrix 
is diagonalised via the bi- unitary transformation M^+ = U*XV\ where U 
and V are complex unitary matrices. The mass matrix for the neutralinos in 
the (B, W, Hi, H 2 ) basis is given by 



Y 



/Mi -M z c p s w M z s p s w \ 

M 2 MzCpCw -M z spc w 

-Mzcpsw M z cpc w -fi 

\ M z sps w -M z spc w -\i ) 



(3) 



where sw(cw) is the sin(cos) of the electro-weak mixing angle 6w Since Y 
is complex symmetric, its diagonalisation requires only one unitary matrix 
N*YN\ 

We study the pair production of charginos, Xi > 



N, via M^o 



(4) 
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Figure 1: Tree- level diagrams for the production of charginos xt anc ^ Xi a ^ 
the LC. 



including effects due to polarised beams. Neglecting the electron-Higgs cou- 
plings, this process occurs at leading order via three diagrams, as seen in 

fig-Hi 



The transition matrix element can be written as 18 



AV(e + e -). xfxj ) = -Q*i3 [v{e + ) lfl P a u{e )] [u(xj )YPpv{xt)] , (5) 



where Q a p denotes the bilinear charges, a = L,R refers to the chirality of the 
e + e _ current and /3 = L, Rto that of the xtxj current. The summation over 
a and (5 is implied. The bilinear charges are comprised of the propagators 
and couplings 



L 



Qll — C-h — — DzGrC~+ ... 
Qrl = C-h — — DzGrC^+~-„. 

iR i n r> ( nR \ * i 1 n I r^R \ nR 



Qlr = C~+ — +D z G L [C H + - 7 ) +— Dp C* - Ci 
e««=^- j7 + Z) z G R (c«_ JZ )*, (6) 

for which the required MSSM couplings for the x^~x~7, Xi'XjZ an d ei> e x/" 
vertices are given by 



ie In 1 



= ( S W$ij ~ U* x Uii - -U* 2 U i2 

C? + -- z =Cf' + „ z (i7->0 5 



C "+-= VJi, (7) 



s w 
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and Gl, Gr, Dz and D„ are defined via 



/ ■> _ s w 2 n — Sw 

Lri — , <-f_R — , 

S W C W C W 

D 7 = — D B 



In the equations above, e denotes the electric charge, m e and are the 
masses of the electron and Z boson. Dz and D„ refer to the propagators of 
the Z boson and sneutrino (of mass mp), in terms of the Mandelstam vari- 
ables s and t. We can neglect the non-zero Z boson width for the considered 
energies. 

One can therefore express the transition matrix element in terms of M 2 , \i 
and tan /3, in addition to the known SM parameters. However, the expected 
accuracy of the measurements at the linear collider is such that one-loop 
corrections become relevant, and we shall see in the following section how 
the higher order expressions depend on many additional MSSM parameters. 



3 NLO contributions and renormalization 



We have calculated the full one-loop corrections to the process e + e — > X1X1 j 
following ref. 19 , within the complex MSSM; examples for the contributing 



self-energy, vertex and box diagrams are shown in fig. [2] For the calculation 
we have used the program FeynArts 20-24 , which allowed an automated 



generation of the Feynman diagrams and amplitudes. Together with the 
packages FormCalc 25-27 and LoopTools 25 we derived the final matrix 



elements and loop integrals. We assume a unit CKM matrix. We regularise 

which ensures that SUSY is preserved, 
In order to obtain finite 



using dimensional reduction 28 30 



via the implementation described in refs. 25 , 31 



results we calculated the counterterms for the involved 2- and 3-point vertices, 
according to the renormalization scheme outlined below, and implemented 
these in the model files of FeynArts. 

A number of one-loop calculations in the gaugino-higgsino sector can be 
found in the literature, mainly in the CP-conserving MSSM §[32}g0], but 
some of which apply a renormalization scheme that is also applicable for 
complex parameters [9 38 . CP-odd observables have also been calculated at 
the one-loop level, for instance in refs. 41 43 , but no dedicated renormal- 



ization scheme was required in these the observables studied were 

UV-finite. Since we intend to extend the current study to the case of com- 
plex parameters, here we follow the approach of refs. 14, 19 closely, where 



a dedicated on-shell renormalization scheme for the chargino and neutralino 
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Figure 2: Examples of one-loop self-energy (upper), vertex (middle) and box 
(lower) diagrams for the production of charginos xt an d Xi a t the LC. 



sector of the MSSM with complex parameters was developed. In the following 
we will therefore only discuss the parameter renormalization of the chargino 
and neutralino sector briefly, and for further details about the chargino field 
renormalization and the renormalization of other sectors we refer the reader 
to refs. 



14,15,19,44 



The mass matrix in the chargino sector, eq. (|2]), is renormalized via 

X ->■ X + 5X, (9) 



where SX is defined by 

sx = 



5M 2 V25{M w sp) 
\f25{M w Cf}) Sfj, 



(10) 



containing the renormalization constants 8M2 and Sfi, as well as renormaliza- 
tion constants (RCs) from other sectors, Sep, 5s p (which can be expressed in 



terms of S tan /?), and 5Mw, defined in ref. 19 . The neutralino mass matrix, 
eq. (|3j) , is similarly renormalized via 



Y ->■ Y + SY, 



5 



where 5Y is defined analogously to SX in eq. (10) and contains the additional 
6Mx, cf. eq. Q. 



Following e.g. ref. 14 , 5Mi, 5M 2 and S/i are determined by choosing 
three out of the total six physical masses of the charginos and neutralinos 
to be on-shell, i.e. the tree-level masses, m^., coincide with the one-loop 
renormalized masses, M%. = + Am^, 

A% = -^Re[^(m|J + E£(m|J] - ^Re[Ef (m|J + Ef (m|J] 

= 0. (12) 



Note that we define the coefficients Y^- R {jp 2 ) and E^ L (p 2 ) of the self energ; 



y 



via 



(p 2 ) = 2/PlSJ.( P 2 ) + ^Eg(p 2 ) + P*EgV) + P R Sf (p 2 ), (13) 

and define the left and right handed vector and scalar coefficients, E^ (p 2 ) 

and flf J L ^ SR (p 2 ) respectively, of the renormalized self-energy analogously. 

As stated earlier, we consider the parameter renormalization as for the 
complex MSSM, such that our setup is easily adaptable for future extensions. 



In ref. 15 19 , it was shown that in the CP violating case, the 1-loop correc- 
tions to the phases of Mi and fi, i.e. <pMi and M respectively are UV finite. 
Therefore we take the approach that these phases can be left unrenormalized. 
We can then determine the necessary conditions to obtain the absolute values 
| Mi |, IM2I and \p,\, depending on which three physical masses are chosen to 
be on-shell. As we have two external charginos, and in order to easily extend 
our setup to the case of X1X2 production, we assume the NCC scheme with 
Xl, Xi and xf on-shell (IJ@|l9|, such that i' = 1 and i" = 1 and 2. Note 
that in choosing the scheme, it is desirable that the on-shell particles should 
contain significant bino, wino and higgsino components, in order that the 
Mi, M2 and \x parameters are accessible. For the above choice, these condi- 
tions are satisfied for all the scenarios defined in sec. [4], in which the lightest 
neutralino always has a sizeable bino-like component. The parameters in 
question of the chargino mass matrix can then be renormalized via, 

8\M X \ = — rjr- ^— ((2Re(e- i ^N i3 N H )Re(U n V n ) 

+ ReiV 2 2 Re {e^U j2 V j2 ))C k + (Re (U n V n )Re (e~^U k2 V k2 ) 
- Re (e-^U j2 V j2 )Re (U kl V kl ))Ni - (ReiV 2 2 Re (e- l ^U k2 V k2 ) 

+ 2Re {e-^"N a N a ) Re (U kl V kl ))C^ , (14) 
1 We adopt the convention that M 2 is real. 



6 



e 



e 



Figure 3: Counterterm diagrams in the MSSM for the production of charginos 
Xi and Xi a t the LC. 

5\M 2 \ =^(Re(e~^U j2 V, 2 )C k -Re(e- i ^U k2 V k2 )C,>), (15) 

5\fi\ = - j (Re (U n V n )C k - Re (U kl V kl )C^ , (16) 
where F, Ci and Ni are defined by 

F =2 (Re (U kl V kl )Re {e^U j2 V j2 ) ~ Re (E^i) Re (e-^U k2 V k2 )) , (17) 

c =Re K+feK+) + e±,«K+)1 + S S*K+) + E SK+)1 

-^25X jfc Re(^y ifc ), (18) 

fe=l,2 

tf< =Re Ko[S^(m|o) + S£,(m|o)] + Egj(ro|o) + E s |(m|o)] 

-^^-fcRe^-JVi*), (19) 

4=1,2 

and the subscripts ± and identify the coefficients of the chargino and neu- 
tralino self-energy respectively. 

Finite results for the process of interest at one-loop are obtained by adding 
the counterterm diagrams shown in fig. [3j Although FeynArts generates 
these diagrams, expressions for the counterterms which renormalize the cou- 
plings defined at tree-level in eq. ^ are required as input, and therefore we 
provide expressions for these here explicitly. For the ■jxt'Xj > ZxfxJ an d 
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eu e xt vertices, these can be expressed as follows, 

5C L + „- =C L + ^ (5Z e + 6 ^)+C L + ^J-^ + -(5Zt l1 +5Z^ ii ), 
xjx, 7 xjx, 7 V 2 / x t x i z 2 2 v , J ' Jj 



x+x, 2 X + X] z y Cw Sw 2 J t 2 



n=l,2 



where the analogous right-handed parts are obtained by the replacement 
L — > R, and 

SC R + „_ = C R + - (V - — + - (SZ Pe + 5Z?) ) 

+ -(c R + -5Z R U + C R + -8ZSJ) . (21) 

Note that the renormalization constants of the SM fields, i.e. Zyy (V = 7, Z) 
and SZg for the vector bosons and electron, and parameters, i.e. SZ e and 
6cw{sw) f° r the electric charge and cos(sin) of the weak mixing angle re- 
spectively, can be found in ref. [19]. The renormalization of the chargino 
fields are defined in the most general way, making use of separate RCs for 
the incoming and outgoing fields, i.e. coefficients SZ± R and SZ± R respec- 
tively for left and right-handed charginos as given in ref. (19). Finally, the 
counterterm for the sneutrino self energy, required for the FeynArts model 
file, takes the form 

/1 777,? \ 

SOb^ = iki l^(6Z 0i + 6Z;y - 5m% - -f(8Z~ n + SZft) , (22) 

for Vi = i) e , u„, i> T , where the sneutrino field and mass RCs, 5Z~. and Sm^-, 



are also defined following ref. 19 



Inital and final state soft radiation must also be included to obtain an 
infra-red finite result as the incoming and outgoing particles are charged, and 



this is done as described in detail in ref. 19 , using the phase-space slicing 
method to define the singular soft and collinear contributions in the regions 
E < AE and 9 < AO respectively. In the soft and collinear limit, regularising 
using electron and photon masses respectively, these can be factorized into 
analytically integrable expressions proportional to the tree-level cross-section 
cr trcc (e + e~ — > XiXi)- However the result is cut-off dependent (i.e. on AE 
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and A9), and removing this dependence requires a calculation of the cross 
section for the three body final state, excluding the soft and collinear regions, 
which can be performed using FeynArts and FormCalc. It is also possible 
to include soft photon radiation in the result for the cross-section obtained 
from FormCalc, however the collinear contribution must be added explicitly, 



for which expressions are given in ref. 19 



4 Fit strategy and numerical results 

4.1 Obtaining MSSM parameters from the fit 

With the loop corrections calculated as in section |3j we can determine the 
fundamental parameters of the MSSM at NLO. In the chargino and neu- 
tralino sectors there are four real parameters, see sec. |2j which we fit to, 

M u M 2 , n, tan/3. (23) 

Additionally, if the scenario is not such that it would already have been 
measured at the LC, we fit to the sneutrino mass, as this enters at tree level 
and will therefore significantly affect cross sections and forward-backward 
asymmteries. At the loop level, the full MSSM parameter set will contribute. 
Depending on the scenario, only limited knowledge about some of these may 
be available. In particular stop sector parameters are likely to be difficult to 
obtain from LHC data and direct production at the LC might not be possible. 
However, our analysis also offers unique sensitivity to these parameters at the 
LC, as stops could significantly contribute to chargino/neutralino observables 
at NLO. 

At the LC, masses are expected to be measured with high precision using 
different methods II]. In the following we use the errors which could be 
achieved using the threshold scan method, however we also investigate how 
the fit precision changes if the masses were obtained from the continuum. In 
case of the cross sections, the experimental uncertainty is dominated by the 



statistical uncertainty 45 



- (24) 
a b 

where 5* and B are the signal and background contributions, respectively. 
In the following, we assume that the statistical uncertainties for the cross 
sections correspond to an integrated luminosity of C = 200 fb _1 per polar- 
isation assuming the efficiency of e = 15%, which includes branching ratios 



9 



for semileptonic final states and selection efficiency of 50% 45 . Similarly, 
for the forward-backward asymmetry we have 



I- A 2 



(25) 



and the total number of events N = N + + iV_ 45 

In order to estimate the theoretical error, we include the uncertainty 
due to the contribution of other unknown MSSM parameters, e.g. the heavy 
pseudoscalar Higgs boson mass m A o, unless explicitly included in the fit. 
The other important source of theoretical uncertainty are NNLO corrections. 
At present, these are only known for chargino and neutralino masses, for 



which the leading SUSY-QCD NNLO correction was calculated in ref. 39 ,40 



Based on these results we estimate the uncertainty due to NNLO corrections 
to be of the order of 0.5 GeV, i.e. comparable to the expected experimental 
uncertainty. Clearly in the future, NNLO results for cross-sections should 
ideally be included in order to match the expected experimental uncertainty. 



We perform a multi-dimensional x fit using Minuit 46 , 47 



X 



O, - €>,, 



(26) 



where the sum runs over the input observables Oi, depending on the scenario, 
with their corresponding experimental uncertainties SO^. 



4.2 Scenarios studied and motivation 



We carry out the fit for three scenarios, SI, S2 and S3, shown in tables [T] 
and [2l Due to the current status of direct LHC searches 48,49 , in all 



Due to the current status of direct LHC searches 
scenarios we require heavy first and second generation squarks and gluinos. 
On the other hand the stop sector is assumed to be light, and in fig. [4} 
for each of these scenarios, the mass corrections for neutralinos x\ an d X3 
are seen to be sensitive to the stop mixing angle. Further, indirect limits 
lead us to choosing mixed gaugino higgsino scenarios, favoured by the relic 



density measurements 50 and relatively high pseudoscalar Higgs masses, in 



light of flavour physics constraints, i.e. the branching ratio B(b — > S7) and 
the anomalous magnetic moment of the muon A(<j^ — 2)/2, all of which we 
checked using mi cr Omegas 2.4.1 



51 52 



Our scenarios are chosen in order to assess the sensitivity to the de- 
sired parameters in a number of possible situations. For example, in S2 we 
study the sensitivity of the fit to large values of M 2 , such that the wino-like 
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Parameter 


Value 


Parameter 


Value 




125 


M 2 


250 


n 


180 


M A o 


1000 


M 3 


700 


tan j3 


10 


M ei>2 


1500 




1500 


M k 


1500 


M qit2 


1500 




400 




650 



Table 1: Table of parameters (with the exception of tan/3 in GeV), for sce- 
narios 1 (SI) and 2 (S2), where S2 is identical to SI except M 2 = 2000 GeV. 
Here M^/^. and M( e / U ) i represent the left and right handed mass parameters 
for of a slepton/squark of generation % respectively, and Af is the trilinear 
coupling for a sfermion /. 



Parameter 


Value 


Parameter 


Value 




106 


M 2 


212 


n 


180 


M A o 


500 


M 3 


1500 


tan j3 


12 


M ei , 2 


125 


M e3 


106 


M h 


180 




1500 


M U3 


450 




-1850 



Table 2: Table of parameters (with the exception of tan f3 in GeV) for scenario 
3 (S3). 



chargino and neutralino are heavy and decoupled from the bino and higgsino- 
like particles. Further, in S1/S2 we consider the case that the sleptons and 
pseudoscalar Higgs bosons are at the TeV scale, and in S3 the case that they 
are relatively light. Note that S3, see tab. |2j is also chosen to be 125 GeV 



such that it is compatible with the recent Higgs results from the LHC 53 ,54 



4.3 Results for scenario 1 

In this scenario, only the charginos and three neutralinos will be accessible 
at the LC. It should be possible to probe the supersymmetric QCD sector, 
with masses of ~1.5 TeV, at the LHC. As input for the fit we therefore use: 
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-1.0 -0.5 0.0 0.5 1.0 

cos#, 




-20: \ : 
-25 r \ /' : 

-1.0 -0.5 0.0 0.5 1.0 

cosd, 

Figure 4: One-loop corrections to the masses of neutralinos \2 (upper) and 
%3 (lower) as a function of the stop mixing angle cos 8t, for scenarios SI 
(blue), S2 (red, dashed) and S3 (green, dotted). 



the masses of the charginos (xf, xf) an d three lightest neutralinos 

(xl,xlxi) 

the light chargino production cross section <j(xtx.i) with polarised 
beams at y/s = 350 and 500 GeV 

the forward-backward asymmetry Aps at = 350 and 500 GeV 



the branching ratio B{b — > sj) calculated using micrOmegas 2.4.1 51 



52 



The input variables, together with errors are listed in tab. |3j Note that 
B{b — > 57) is included in order to increase sensitivity to the third generation 
squark sector, and the estimated experimental precision of 0.3 • 10^ 4 , taken 
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Observable 



Tree value Loop corr. Error exp. Error th. 



m xf 


149.0 




U.l (U.ZJ 






292.3 




0.5 (2) 






106.9 




0.2 




X2 


164.0 


2.0 


0.5 (1) 


0.5 


X3 


188.6 


-1.5 


0.5 (1) 


0.5 


a (XlXl )(-0.8,0.6) 


2347.5 


-291.3 


8.7 


2.0 


a (Xl Xl )(0.8 -0.6) 


224.4 


7.6 


2.7 


0.5 


a (XlXl ) (-0.8,0.6) 


1450.6 


-24.4 


6.7 


2.0 


a {Xl Xl ) ^0.8,-0.6) 


154.8 


12.7 


2.0 


0.5 


Af° B (%) 


-2.2 


6.8 


0.8 


0.1 


Af° B (%) 


-2.6 


5.3 


1.0 


0.1 



Table 3: Observables (masses in GeV, cross sections in fb) used as input for 
the fit in SI, tree- level values and loop corrections are specified. Here the 
superscript on a and A FB denotes in GeV, and the subscript on a denotes 
the beam polarisation (V(e~), V(e + )). The central value of the theoretical 
prediction, B{b — > 57) = 3.3 • 10~ 4 GeV, calculated using state-of-the-art 



tools 51,521, is also included in the fit. Errors in brackets are for masses 



obtained from the continuum. See text for details of error estimation. 



from ref. 55 , is adopted. We found that the impact of the muon anomalous 
magnetic moment is negligible in this scenario, mainly due to the heavy 
smuon sector. 

In SI we fit 8 MSSM parameters: Mi, M 2 , /i, tan/3, m^, cos 6^-, , and 
mi 2 . We explicitly check that the impact of other parameters, eg. the masses 
of the 1st and 2nd generation squarks, is negligible. The small dependence 
on the A mass is included as an additional source of error. Note that there 
are no theoretical errors for masses chosen to be on-shell. These masses can 
be related to the fundamental parameters exactly via tree level relations, and 
are included in the fit. 

The results of the fit are given in tab. |4} We find that the gaugino and 
higgsino mass parameters are determined with an accuracy better than 1%, 
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Parameter Threshold fit Continuum fit 



M 1 


125±0.3 


(±0.7) 


125±0.6 


(±1.2) 


M 2 


250±0.6 


(±1.3) 


250±1.6 


(±3) 


/' 


180±0.4 


(±0.8) 


180±0.7 


(±1.3) 


tan j3 


10±0.5 


(±1) 


10±1.3 


(±2.6) 


mc 


1500±24 




1500±20 


(±40) 


cos 61 


0±0.15 


( +0A ) 






m- h 




/at limit \ 
vat limit/ 






m h 


800±?~ 


/+1000\ 
1-290 ) 


800^5° 


/at limit \ 
Vat limit/ 



Table 4: Fit results (masses in GeV) for SI, for masses obtained from thresh- 
old scans (threshold fit) and from the continuum (continuum fit). Numbers 
in brackets denote 2a errors. 



while tan /3 is determined with an accuracy of 5%. Excellent precision of 2-3% 
is obtained for the otherwise unobservable sneutrino. Including NLO effects 
even allows us to constrain the parameters of the stop sector. Although the 
precision shown in tab. [4] is rather limited, this could lead to an important 
hint concerning the masses of the stops, allowing for a well-targeted search 



at the LHC This could be another example of LC-LHC interplay 56 . 

Finally, in tab. [4] we compare the fit results using masses of the charginos 
and neutralinos obtained from threshold scans and from the continuum. For 
the latter, the fit quality is seen to deteriorate, with errors on the fundamental 
parameters almost doubling, clearly indicating the need to measure chargino 
and neutralino masses via threshold scans. 



4.4 Results for scenario 2 

In this scenario, where the M 2 parameter is set to 2 TeV, only the light 
chargino and three lightest neutralinos will be accessible at the LC. As input 
for the fit we therefore use: 

• the masses of the lighter chargino (x%) an d neutralinos {X11X21X3) 

• the light chargino production cross section <j(xtXi) with polarised 
beams at y/s = 400 and 500 GeV 
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0.0 
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1.0 
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^(%) 


0.0 


5.0 
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0.1 



Table 5: Observables (masses in GeV, cross sections in fb) used as an input 
for the fit in S2, as in tab. |5j The central value of the theoretical prediction, 



B(b — > sj) = 3.3 ■ 10 4 GeV, calculated using state-of-the-art tools 51 52 
is also included in the fit. See text for details of error estimation. 



• the forward-backward asymmetry Afb at = 400 and 500 GeV 

• the branching ratio B(b — > 57). 

As we again find that the muon anomalous magnetic moment has a negligible 
impact, it is not used in the fit. The input variables, together with errors are 
listed in tab. [5l 

We again fit 8 MSSM parameters: Mi, M 2 , /1, tan/3, ma, cos Of, m^, and 
?7i£ 2 . The impact of other parameters, except the heavy Higgs boson mass, 
can be neglected. The results from the fit are given in tab. |6j The higgsino 
and bino mass parameters are well constrained in this scenario since bino-like 
neutralino and all higgsinos are directly accessible. Even though the winos 
are not directly accessible, the wino mass parameter M2 can be constrained 
with 10% accuracy at la level. An accuracy of 20% is achieved for tan/3, 
significantly worse than in SI. This can be understood by the fact that the 
mixing in S2 between chargino states is weak due to M 2 being heavy, and the 
constraint on tan f3 is dependent on this mixing. No limits can be derived on 
the sneutrino mass, due to the Yukawa suppressed coupling of the higgsino- 
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Parameter Fit result 
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m h 


n +0.13 /+0.4\ 
u -0.09 V-0.3J 

8nn+ 300 i'+soon 

ouu -200 1-400J 



Table 6: Fit results (with the exception of tan j3 and cos^ in GeV) for S2, 
as in tab. |4[ where numbers in brackets denote 2a errors. 



like xt to the electron and sneutrino. We are however, as shown in tab. |6j 
still able to derive limits on the stop masses and mixing parameter. 

4.5 Results for parameters in scenario 3 

This final scenario features the richest phenomenology of the studied bench- 
mark scenarios. As input for the fit we therefore use: 

• the masses of the charginos X2 ) 3Jid neutralinos (XijX^Xs) 

• the light chargino production cross section cr(xfXi) with polarised 
beams at y/s = 400 and 500 GeV 

• the forward-backward asymmetry Afb at = 400 and 500 GeV 

• the Higgs boson mass, rrih 

• the branching ratio B(b — > 57) 

• the anomalous muon magnetic moment 

Compared to the previous scenarios, these observables are supplemented by 



the Higgs boson mass, m h , calculated using FeynHiggs 2.9.1 57 60]. The 
estimated experimental precision at the LC for m^, taken from ref. fl], is 
adopted. We further assume the theoretical uncertainty on the Higgs boson 
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Table 7: Observables (masses in GeV, cross sections in fb) used as an input 
for the fit in S3, as in tab. [3] The central theoretical predictions B(b — > 
si) = 2.7 • 10" 4 , A(g„ - 2)/2 = 2.4 • 10" 9 and m h = 125 GeV, calculated 
using state-of-the-art tools (51 52 57 60 , are also included in the fit. See 



text for details of error estimation. 



mass to be 1 GeV [60]. As before, the remaining two observables, the 
branching ratio B(b — > 57) and the anomalous muon magnetic moment are 



calculated using micrOmegas 2.4.1 51,52 , and a projected experimental 
error on the anomalous muon magnetic moment of 3.4- 10 -10 is employed 61 



which we assume would dominate over the theoretical uncertainty. The input 
variables, along with errors, are summarised in tab.[7| Because the sneutrino 
is now directly accessible, we assume that its mass has been measured. 

In addition to the 7 MSSM parameters fitted in scenarios 1 and 2: Mi, 
Mi, \x, tan (3, cos Of, m^, and m 4 ~ 2 , we now explicitly include the heavy Higgs 
boson mass in the fit and remove the sneutrino mass. The results of the fit 
are collected in tab. [8] The parameters of the electroweak gaugino-higgsino 
sector are determined with high precision. Due to a significant mixing in the 
stop sector, and the improvement in the fit quality due to the inclusion of the 
higgs mass, we find that the fit is now also sensitive to the mass of the heavy, 
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Parameter Fit result 
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Table 8: Fit results (with the exception of tan/3 and cos^ in GeV) for S3, 
as in tab. |4| but in addition including results for the masses of the heavier 
stop mass (m^ ) and the pseudoscalar higgs boson (m A o). 

left-handed stop. The accuracy is better than 20% for this particle even 
though it is far beyond the reach of the LC and also LHC In addition, an 
upper limit on the mass of the heavy Higgs boson can be placed at 1000 GeV, 
at the 2cr level. 

5 Conclusion 

The evidence for the Higgs boson and dark matter, when examined in the 
context of supersymmetry, suggest the possibility of a light /i and Mi. We 
have extended previous analyses, which fitted observables for chargino pro- 
duction at the LC to extract fundamental MSSM parameters, by incorporat- 
ing NLO corrections. The loop corrections are calculated for all observables 
fitted, namely the polarised cross-sections and forward backward asymmetry 
for chargino production as well as the X11X2 an d Xi,X2,X3 masses, in an 
on-shell scheme which facilitates the extension to the complex case. We have 
fitted these observables for three complementary scenarios. We found that 
on including NLO corrections, when M 1; M 2 and /1 are light they can be 
determined to percent-level accuracy, and tan (3 to < 5%. Further we showed 
that obtaining masses of the charginos and neutralinos from the continuum 
as opposed to threshold scans would result in the uncertainty on the funda- 
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mental parameters almost doubling, reinforcing the importance of threshold 
scans for mass measurements. As a heavy M 2 is still a viable possibility, 
we also considered M 2 = 2000 GeV, and found that the sensitivity to M 2 
is approximately 10%. As the error on tan j3 is dependent on the degree of 
mixing in the chargino sector, here it increases to ~ 20%. Note that the 
inclusion of B(b — » 57), as well as the use of masses determined via threshold 
scanning, in the fit was seen to improve the sensitivity to the stop sector. We 
finally considered a scenario compatible with the latest Higgs results. For 
this scenario we found that including B{b — > 57), A(p At — 2)/2 and m h in the 
fit, along with the significant mixing in the stop sector, helped to obtain an 
accuracy better than 20% on the mass of the heavy, left-handed stop, even 
though this particle is far beyond the reach of the LC and also the LHC We 
also included m A o in the fit, and found that it would be possible to place a 
2a upper bound on this parameter of 1000 GeV. In summary, we have shown 
that incorporating NLO corrections is required for the precise determina- 
tion of the fundamental parameters of the chargino and neutralino sector at 
the LC, and could further provide sensitivity to the parameters describing 
particles which contribute via loop corrections. 



Acknowledgements 

The authors gratefully acknowledge support of the DFG through the grant 
SFB 676, "Particles, Strings, and the Early Universe" , as well as the Helmholtz 
Alliance, "Physics at the Terascale". This work was also partially sup- 
ported by the Polish National Science Centre under research grant DEC- 
2011/01/M/ST2/02466 and the MICINN, Spain, under contract FPA2010- 
17747; Consolider-Ingenio CPAN CSD2007- 00042. KR thanks as well the 
Comunidad de Madrid through Proyecto HEPHACOS S2009/ESP-1473 and 
the European Commission under contract PITN-GA-2009-237920. 

References 

[1] ECFA/DESY LC Physics Working Group Collaboration, 
J. Aguilar-Saavedra et al. |hep- ph/0 106315, 

[2] ILC Collaboration Collaboration, E. Brau, James et al. 
larXiv: 0712 .19501 



[3] ILC Collaboration, G. Aarons et al. 



arXiv: 0709. 1893 



[4] H. Goldberg Phys. Rev. Lett. 50 (1983) 1419. 



19 



[5] J. R. Ellis, J. Hagelin, D. V. Nanopoulos, K. A. Olive, and 
M. Srednicki Nucl.Phys. B238 (1984) 453-476. 

[6] L. J. Hall , D. Pinne r, and J. T. Ruderman JEEP 1204 (2012) 131, 
[arXiv: 1112 .2703] . 

[7] F. Brummer and W. Buchmuller JEEP 1107 (2011) 010, 
[arXiv: 1105 .0802] . 

[8] K. Desch, J. Kalinowski, G. A. Moortgat-Pick, M. Nojiri, and 
G. Polesello JEEP 0402 (2004) 035, |hep -ph/ 03 12069] . 

[9] W. Oiler, H. Eb erl, and W. Majerotto Phys.Rev. D71 (2005) 115002, 
[hep-ph/0504109| . 

[10] T. Fritzsche and W . Hollik Nucl.Phys. Proc.Suppl. 135 (2004) 102-106, 
[hep-ph/0407095| . 

[11] W. Kilian, J. Reu ter, and T. Robens Eur.Phys.J. C48 (2006) 389-400, 
[hep-ph/0607127 . 

[12] T. Robens, J. Kalinowski, K. Rolbiecki, W. Kilian, and J. Reuter Acta 
Phys.Polon. B39 (2008) 1705-1714, jarXiv: 0803. 4161 . 

[13] T. Fritzsche. PhD thesis, Universitaet Karlsruhe, 2005. 

[14] A. Fowler and G. W eiglein JEEP 1001 (2010) 108, 
[arXiv: 0909 .5165] . 

[15] A. Fowler. PhD thesis, Durham University, 2010. 

[16] A. Chatterjee, M. Drees, S. Kulkarni, and Q. Xu |arXiv: 1107.5218 

[17] H. E. Haber and G. L. Kane Phys.Rept. 117 (1985) 75-263. 

[18] S. Choi, A. Djouadi, M. Guchait, J. Kalinowski, H. Song, et al. 
Eur.Phys.J. C14 (2000) 535-546, lhep-ph/0002033] . 

[19] A. Bharucha, A. Fowler, G. Moortgat-Pick, and G. Weiglein 
larXiv: 1211 .31341 

[20] J. Kublbeck, M. Bohm, and A. Denner Comput.Phys.Commun. 60 
(1990) 165-180. 

[21] A. Denner, H. Eck, O. Hahn, and J. Kublbeck Nucl.Phys. B387 (1992) 
467-484. 



20 



[22] J. Kublbeck, H. Eck, and R. Mertig Nucl.Phys.Proc.Suppl. 29A (1992) 
204-208. 

[23] T. Hahn Comput. Phys. Commun. 140 (2001) 418-431, 
|hep-ph/0012260| . 

[24] T. Hahn and C. Schappacher Comput.Phys. Commun. 143 (2002) 
54-68, [ hep-ph/0105349| . 

[25] T. Hahn and M. Perez- Victoria Comput.Phys. Commun. 118 (1999) 
153-165, [hep-ph/9807565] . 

[26] T. Hahn Comput. Phys. Commun. 178 (2008) 217-221, 
[hep-ph/0611273| . 

[27] T. Hahn and M. R auch Nucl.Phys.Proc.Suppl. 157 (2006) 236-240, 
|hep-ph/0601248| . 

[28] W. Siegel Phys.Lett. B84 (1979) 193. 

[29] W. Siegel Phys.Lett. B94 (1980) 37. 

[30] D. Stockinger JEEP 0503 (2005) 076, |hep-ph/0503129] . 

[31] F. del Aguila, A. Culatti, R. Munoz Tapia, and M. Perez- Victoria 
Nucl.Phys. B537 (1999) 561-585, |hep-ph/980645l] . 

[32] A. Lahan as, K. Tamvakis, a nd N. Tracas Phys.Lett. B324 (1994) 
387-396, |hep-ph/93 12251] . 

[33] D. Pierce and A. P apadopoulos Phys.Rev. D50 (1994) 565-570, 
|hep-p h/9312248| . 

[34] D. Pierce and A . Papadopoulos Nucl.Phys. B430 (1994) 278-294, 
|hep-ph/9403240| . 

[35] H. Eberl, M. K incel, W. Majerott o, and Y. Yamada Phys.Rev. D64 
(2001) 115013, |hep-ph/0104109] . 

[36] T. Fritzsche and W . Hollik Eur.Phys.J. C24 (2002) 619-629, 
|hep-ph/0203159| . 

[37] W. Oiler, H. Eb erl, W. Majerott o, and C. Weber Eur.Phys.J. C29 
(2003) 563-572, [hep-ph/0304006] . 



21 



[38] M. Dre es, W. H ollik, and Q. Xu JEEP 0702 (2007) 032, 
|hep-ph/0610267| . 

[39] R. Schofbeck and H. Eberl Phys.Lett. B649 (2007) 67-72, 
|hep-ph/0612276| . 

[40] R. Schofbeck and H . Eberl Eur.Phys.J. C53 (2008) 621-626, 
jarXiv: 0706 .0781] . 

[41] K. Rolbiecki and J Kalinowski Phys.Rev. D76 (2007) 115006, 
[arXiv: 0709. 2994] . 

[42] H. Eberl, T. Gajdosik, W. Majerotto, and B. Schrausser Phys.Lett. 
B618 (2005) 171-181, [hep-ph/050"2112] . 

[43] P. Osland and A. Vereshagin Phys.Rev. D76 (2007) 036001, 



arXiv: 0704. 2165 



[44] A. Bharucha, S. Heinemeyer, F. von der Pahlen, and C. Schappacher 
larXiv: 1208 .41061 

[45] K. Desch, J. Kalinowski, G. Moortgat-Pick, K. Rolbiecki, and 
W. Stirling JEEP 0612 (2006) 007, |hep-ph/0607104] . 

[46] F. James and M. Roos Comput.Phys.Commun. 10 (1975) 343-367. 

[47] F. James CERN Program Library Long Writeup D506 (1994). 

[48] ATLAS Collaboration Collaboration, G. Aad et al. 
larXiv: 1208.09491 

[49] CMS Collaboration Collaboration, S. Chatrchyan et al. 
larXiv: 1207. 17981 

[50] WMAP Collaboration Collaboration, E. Komatsu et al. 
Astrophys.J.Suppl. 192 (2011) 18, [arXiv: 1001.4538| . 

[51] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov 



Comput.Phys.Commun. 176 (2007) 367-382, |hep-ph/0607059 



[52] G. Belanger, F. Boudjema, P. Brun, A. Pukhov, S. Rosier-Lees, et al. 
Comput.Phys.Commun. 182 (2011) 842-856, |arXiv: 1004. 1092| . 

[53] ATLAS Collaboration Collaboration, G. Aad et al. Phys.Lett. B716 
(2012) 1-29, |arXiv:1207.7214|. 



22 



[54] CMS Collaboration Collaboration, S. Chatrchyan et al. Phys.Lett. 
B716 (2012) 30-61, jarXiv: 1207.7235] . 

[55] SuperB Collaboration Collaboration, B. O'Leary et al. 
larXiv: 1008 .15411 

[56] LHC/LC Study Group Collaboration, G. Weiglein et al. Phys.Rept. 
426 (2006) 47-358, lhep-ph/04 10364] . 

[57] S. Heinemeyer, W. Hollik, and G. Weiglein Comput.Phys.Commun. 
124 (2000) 76-89, [hep-ph/9812320] . 

[58] S. Heine meyer, W. Hollik, and G. Weiglein Eur.Phys.J. C9 (1999) 
343-366, [hep -ph/9812472] . 

[59] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich, and G. Weiglein 
Eur.Phys.J. C28 (2003) 133-143, |hep-ph/02 12020 . 

[60] M. Frank, T. Halm, S. Heinemeyer, W. Hollik, H. Rzehak, et al. JHEP 
0702 (2007) 047, |hep-ph 706Tll26] . 

[61] R. Carey, K. Lynch, J. Miller, B. Roberts, W. Morse, et al. 



23 



